Existence of an upper limit on the density of excitons in carbon nanotubes by 
diffusion-limited exciton-exciton annihilation: Experiment and theory 
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Through an investigation of photoemission properties of highly-photoexcited single- waUed carbon 
nanotubes, we demonstrate that there is an upper limit on the achievable excitonic density. As the 
intensity of optical excitation increases, all photoluminescence emission peaks arising from different 
chirality single-walled carbon nanotubes showed clear saturation in intensity. Each peak exhibited a 
saturation value that was independent of the excitation wavelength, indicating that there is an upper 
limit on the excitonic density for each nanotube species. We propose that this saturation behavior 
is a result of efficient exciton-exciton annihilation through which excitons decay non-radiatively. In 
order to explain the experimental results and obtain excitonic densities in the saturation regime, 
we have developed a model, taking into account the generation, diffusion-limited exciton-exciton 
annihilation, and spontaneous decays of one-dimensional excitons. Using the model, we were able to 
reproduce the experimentally obtained saturation curves under certain approximations, from which 
the excitonic densities were estimated. The validity of the model was confirmed through comparison 
with Monte Carlo simulations. Finally, we show that the conventional rate equation for exciton- 
exciton annihilation without taking into account exciton diffusion fails to fit the experimentally 
observed saturation behavior, especially at high excitonic densities. 

PACS numbers: 78.67.Ch,71.35.-y,78.55.-m 
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I. INTRODUCTION 

The optical and electronic properties of low- 
dimensional materials have been an important subject 
of study in the field of condensed matter physics. In 
particular, one-dimensional (1-D) materials are predicted 
to possess unique properties that are distinctly differ- 
ent from those at higher dimensions fii^ primarily due 
to the enhanced Coulomb interactions among the quan- 
tum confined charge carriers. One common feature 
of optically-excited low-dimensional systems is the for- 
mation of strongly bound electron-hole {e-h) pairs, or 
excitons,^ which dominate interband optical spectra. 1- 
D semiconductors are expected to show an almost com- 
plete suppression of optical absorption at the band edges, 
with a significant fraction of the total oscillator strength 
taken by the lowest excitonic state<^i^ 

Early reports of lasing from semiconductor quantum 
wires (QWRs)^.!'^ invoked much interest in the physics 
of high density 1-D excitons. A number of studies have 
thus far been performed on such QWR systems during 
the last two decades to understand many-body phenom- 
ena (e.g., lasing, band-gap renormalization, biexciton for- 
mation, and the Mott transition), ^.>^iiSiiiiiSiiliiiii^ but 
many aspects are still under debate and not well un- 
derstood. More recently, studies on high-density e-h 
pairs have been extended to novel 1-D materials such 
as conjugated polymer a^^'^^i^^'^^ and single- walled car- 
bon nanotubes (SWNTs),^^!^!^^ The latter are tubular 
materials made of sp^-bonded carbon atoms, ^'^ attracting 
much recent interest from diverse research fields due to 
their unique properties. Semiconducting SWNTs are 
known to have extremely strong quantum confinement 



of ~ 1 nm, giving rise to large exciton binding energies 
on the order of 0.5-1 eVr^i^'21 much larger than those 
of GaAs QWRs (~ 20 meV)^'- and comparable to or 
larger than those of conjugated polymers (~ 0.4 eV— 
and < 0.1 eV^). 



Here, we report results of experimental and theoreti- 
cal investigations on the properties of photoluminescence 
(PL) from excitons in SWNTs through nonlinear pho- 
toluminescence excitation (PLE) spectroscopy using in- 
tense optical pulses. From the clear saturation behav- 
ior observed in the intensities of all the PL features as 
a function of excitation laser intensity as well as the 
complete flattening of the PLE spectra observed at very 
high laser intensities, we show the existence of an upper 
limit on the density of excitons that can be accommo- 
dated in SWNTs. Such an upper limit is considered to 
be caused by the diffusive motion of the excitons^'^'^ 
combined with highly rapid and efficient exciton-exciton 
annihilation (EEA) in SWNTs.— As described in Section 
mil we have developed a theoretical model for describing 
diffusion- limited EEA processes in 1-D, which enabled 
us to simulate the PL saturation curves and estimate the 
densities of excitons in SWNTs as a function of excitation 
intensity. 



A portion of this work was described in our earlier let- 
ter.3? The purpose of the present paper is to provide a 
complete description of both experimental and theoreti- 
cal aspects of this study. 
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II. EXPERIMENT 
A. Experimental methods 

The sample was prepared by ultrasonicating CoMo- 
CAT SWNTs in D2O with 1 wt% sodium chelate for 1 
hour, followed by ultracentrifugation at 111,000 g for 4 
hours. This centrifugation condition is sufficient to re- 
move SWNT bundles effectively^^ Only the upper 50 % 
of the supernatant was collected and used for the ex- 
periment. The solution was put in a 1-mm-thick quartz 
cuvette. The optical density of the sample around the 
E22 resonance was below 0.2, which helped avoid non- 
uniform excitation and re-absorption of the emitted PL 
within the sample. The excitation source was an optical 
parametric amplifier (OPA), producing ^ 250 fs pulses 
at a repetition rate of 1 kHz, tunable in the visible and 
near-infrared ranges, pumped by a chirped-pulse ampli- 
fier (Clark-MXR, Inc., CPA-2010). Optical filters were 
carefully selected and set in the beam path to thoroughly 
eliminate any parasitic wavelength components (mostly 
in the ultraviolet and near-infrared regions) contained in 
the OPA beam. The OPA beam was focused onto the 
sample to a spot size of 300-400 fj,m. Only the central 
~ 2 mm portion of the OPA beam profile (~ 6 mm) 
was taken out by using an aperture just before the focus 
to enhance the spot uniformity at the sample. The PL 
from the sample was focused onto the monochromator 
entrance and recorded with a liquid-nitrogen-cooled In- 
GaAs 1-D array detector. The obtained PL spectra were 
corrected for the wavelength-dependence of the grating 
efficiency and detector sensitivity. 

For the data shown in Fig. [TJc), a different sample as 
well as a different excitation light source was used for ver- 
ifying the universality of the phenomena observed. In this 
case, the sample was a dried film of CoMoCAT SWNTs 
embedded in t-carrageenan, formed by drying a mixed gel 
of i-carrageenan and the centrifuged supernatant of Co- 
MoCAT SWNTs on the surface of a sapphire substrate. 
The sapphire substrate served as a mechanical support as 
well as a heat sink of the film during the measurements. 
The excitation light was 1 kHz and ^ 250 fs optical pulses 
with a central wavelength of 653 nm (FWHM = 10 nm), 
produced by filtering whitelight pulses generated by fo- 
cusing the CPA beam onto a sapphire crystal. 

B. Experimental results 

Figure [Ija) compares two PL spectra. The black solid 
curve was obtained using the OPA with a wavelength of 
654 nm (or 1.90 eV) and a pulse energy of 29 nJ, while 
the red dotted curve was obtained using a weak (100 /iW) 
CW laser with a wavelength of 658 nm (or 1.88 eV). It is 
seen that the relative intensities of different PL peaks are 
drastically different between the two curves. The inset 
confirms that the two spectra coincide accurately when 
the OPA pulse energy was kept very low (300 pJ). Fig- 



OPA, 654 nm CW laser diode, 658 nm 




950 1000 1050 1100 1150 
Wavelength (nm) 




950 1000 1050 1100 1150 
Wavelength (nm) 



FIG. 1: (color online) Pump-intensity-dependent PL spectra 
measured for (a, b) the centrifuged supernatant of CoMoCAT 
SWNTs dispersed in D2O and (c) CoMoCAT SWNTs em- 
bedded in a dried t-carrageenan film: (a) Black solid — spec- 
trum obtained with OPA pulse (654 nm, 29 nJ). Red dotted 
curve — spectrum obtained with a CW laser diode (658 nm, 
100 /iW). Inset shows that the two spectra coincide when the 
OPA pulse energy is very low (300 pJ). (b) Change of PL spec- 
tra with the pulse energy of OPA beam (654 nm) varied be- 
tween 1 nJ and 30 nJ in the order of #1 to #7. Curve #4 cor- 
responds to the highest fluence (~L3 x 10^'* photons/cm^), 
(c) Change of PL spectra with the pulse energy of 653 nm 
light (FWHM = 10 nm) varied between 1 nJ and 20 nJ in the 
order of #1 to #7. 



ure [TJb) shows PL spectra measured with pulse energies 
of 1 nJ (curves #1 and #7), 4 nJ (#2 and #6), 10 nJ (#3 
and #5), and 30 nJ (#4). The (7,5) peak is dominant at 
low fiuences while the (6,5) peak becomes dominant at 
high fiuences. It is important to note that the different 
curves were taken in the order of #1 to #7, demonstrat- 
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ing that the observed changes are reproducible and are 
not caused by any laser-induced permanent change in the 
sample. Additionally, note that the PL intensities tend to 
saturate at high laser fluences, while their peak positions 
do not change at all. 

Figure [He) shows PL spectra measured for the dried 
t-carrageenan film using 653 nm optical pulses with a 
FWHM bandwidth of 10 nm at different pulse energies. 
The pulse energies were 1 nJ (curves #1 and #7), 3 nJ 
(#2 and #6), 10 nJ (#3 and #5), and 20 nJ (#4) mea- 
sured in the order of #1 to #7. Figure [TJc) exhibits the 
same behavior as that shown in Fig.[TJb), demonstrating 
that the observed changes shown in Fig. [Ijb) did not re- 
sult from any artifacts, e.g., caused by the fluidic nature 
of the sample or by the unnoticed parasitic wavelength 
components in the OPA beam. In the following, we use 
the excitation pulse fluence in terms of the number of in- 
cident photons per cm^ per pulse to express the intensity 
of excitation pulses. 

Figures [2l^a)-[2jd) show PLE maps taken with vari- 
ous pump fluences. The step size for the pump pho- 
ton energy was 20 meV. The data taken with the low- 
est fluence (1.2 x 10^^ photons/cm^) [Fig. El^a)] is es- 
sentially the same as that taken with low-intensity CW 
light. However, as the fluence is increased [Figs. [H^b)- 
[UJd)], the E22 excitation peaks gradually broaden and 
eventually become completely flat at the highest fluence 
(1.2 X 10^^ photons/cm^) — i.e., PL intensities become 
independent of the excitation wavelength. The corre- 
sponding PLE spectra are shown in Figs. [2l[e)-[2l^h) for 
three PL wavelengths at 983, 1034, and 1125 nm. Again, 
such changes in the PLE spectra were reproducible over 
the fluence range tested here, indicating that no sample 
damage was induced. 

In order to obtain PL intensity (/pl) versus pump in- 
tensity (/pump) relationships for different emission peaks, 
we measured PL spectra at different photon ffuences for 
various excitation wavelengths. Each PL spectrum was 
decomposed and fltted by multiple peaks corresponding 
to the SWNT types/chiralities involved in the measured 
wavelength range. 50% Gaussian -I- 50% Lorentzian line- 
shape was assumed, and the decomposition was per- 
formed by optimizing the peak-width so that the decom- 
position gives the best fltting for the original PL spec- 
trum. The optimum widths at the highest fluence was 
larger by ^ 15 % than those at the lowest fluence in 
Fig. [1] and such an increment of the width is consid- 
ered to be caused by the enhanced interactions among 
excitons or their reduced lifetime in the presence of high 
density excitons. Throughout the decomposition analy- 
sis performed, peak positions of all the PL features and 
the ratios among their widths were flxed regardless of the 
excitation wavelength and fluence. 

Figure [3] shows the obtained integrated PL peak inten- 
sities (/pl) plotted against the incident photon fluence 
(/pump) for (6,5), (7,5), and (8,3) SWNTs at excitation 
wavelengths of 570, 615, and 658 nm. The resonance 
wavelengths of these SWNT types at the £^22 levels are 
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FIG. 2: (color online) Evolution of PLE data with increasing 
pump pulse fluence: (a) 1.2 x 10^^, (b) 1.2 x 10^^, (c) 4.1 x 
10^^, and (d) 1.2 x lO" photons/cm^, (e-h): PLE spectra 
corresponding to (a)-(d) at emission wavelengths of 983 nm 
(circles), 1034 nm (squares), and 1125 nm (triangles). 



approximately 570, 647, and 673 nm, respectively. It can 
be seen that the integrated PL intensity begins to satu- 
rate at a lower (higher) fluence when SWNTs are excited 
resonantly (non-resonantly) . Unexpectedly fast satura- 
tion of the PL from (7,5) with 570 nm excitation (which 
is non-resonant) can be attributed to its proximity to the 
phonon sideband at ~ 585 nmii'* 



C. Data interpretation 

We interpret these observations as results of very ef- 
ficient exciton-exciton annihilation, ^ ^■^^i-'^^'^^i^^i^^ a non- 
radiative process that occurs at high exciton densities 
where two excitons are spatially close enough to interact 
with each other, resulting in the annihilation of the two 
excitons and simultaneous creation of an e-h pair in a 
higher energy state (either as a bound exciton or an un- 
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FIG. 3: (color online) Integrated PL intensity versus pump 
fluence for (6,5), (7,5), and (8,3) SWNTs. Pump wavelengths 
were 570 nm (circles), 615 nm (squares), and 654 nm (trian- 
gles). The error bars account for ±5%. The solid and dashed 
curves are fitting by Eq. (|15fl and Eq. (|20p . respectively. 



bound free e-h pair). We assume that the formation of 
Ell excitons occurs in a very short time scale after an op- 
tical excitation around E22, because of much faster £'22- 
to-i?ii relaxation (e.g., ~40 fs2£) than the duration of our 
OPA pulse (^--^250 fs). Thus, excitons quickly accumulate 
in the En state during and right after photo-creation of 
e-h pairs. However, the number of excitons that can be 
accommodated in the £'11 state is limited by EEA. As 
the exciton density, rix, approaches its maximum value, 
EEA begins to prevent a further increase by efficiently re- 
moving excitons non-radiatively, which explains the PL 
saturation behavior. Since EEA serves as a bottleneck 
for the exciton density, the PL intensity becomes insen- 
sitive to whether the excitons were created resonantly or 
non-resonantly around the E22 level and independent of 



the pump wavelength, resulting in the flattening of PLE 
spectra. Namely, at very high pump fluences, the PL in- 
tensity is determined not by how efficiently excitons are 
created but by how many En excitons can be accommo- 
dated within a particular type of SWNT as well as by the 
relative abundance of that type of SWNT in the sample. 
We also performed optical transmission spectroscopy in 
the E22 range using OPA pulses and found that the ab- 
sorption spectra do not exhibit any change even at high 
pulse fluences. Thus, nonlinear optical effects such as 
phase-space filling in the E22 range are not playing any 
role in the observed PLE broadening/flattening and PL 
saturation. 



III. THEORY 

In this section, we develop a theoretical model for ex- 
plaining the experimental results, taking into account the 
generation, diffusion-limited EEA, and spontaneous de- 
cays of 1-D excitons. Under certain approximations, the 
model provides a direct analytical relationship between 
the intensity of the excitation light (/pump) and that of 
the emitted PL (/pl) for limiting cases, which allows us 
to estimate the density of excitons in SWNTs through fit- 
ting to the experimentally obtained /pump vs. /pL curves. 



A. Model 

Figure Ufa) shows a schematic energy diagram of the 
excitons in semiconducting SWNTs under consideration. 
We are interested in calculating the population of exci- 
tons N (indicated by the dotted box) in the lowest en- 
ergy state £^11. First, excitons are created by optical 
excitation at an energy around the E22 level, which is 
typically in the visible wavelength range. The excita- 
tion intensity is denoted by /pump- As soon as excitons 
are created, they decay to the En level within a very 
short time (~40 fs)— by transferring their energies to 
the lattice. Recent studies have reported that the ex- 
citons created around E22 levels primarily decay to the 
Ell level with a probability close to unityj ^^i'^^ The infiux 
of excitons to the £^11 level is denoted by Gin. On the 
other hand, the spontaneous decay time rtot of the En 
excitons to the ground state (G. S.) has been reported to 
be 10-100 ps.'^^''*'''*! Such a spontaneous decay consists 
of radiative and non-radiative processes, with respective 
rates 7^ and 7nr (s~^), where 7r -I- 7nr = 7tot = 7"tot""^- 
The outflux of excitons from the En level via the spon- 
taneous decay process is denoted by Gout- Therefore, the 
flux of the PL photons or the PL intensity (/pl) is equal 
to ?/Gout, where 77 (= 7r/7tot) denotes the branching ratio 
for the radiative decay from the En level. 

As the density of excitons increases, the EEA process 
becomes important. If the e-h pair created in the higher 
energy state returns back to the En exciton level with 
a probability of A (0 < A < 1), the initial two excitons 
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FIG. 4: (color online) (a) Schematic energy diagram of the 
system considered in the model. The dotted box enclosing 
the lowest energy level (-Eii) is the domain of interest where 
excitons are populated. All the symbols are defined in the 
text, (b) Schematic description of A'' = 4 excitons randomly 
distributed over a SWNT with a length of L. The horizontal 
arrow with a length of denotes the average length traversed 
by one exciton during its spontaneous decay lifetime not . The 
lower part is the equivalent of the upper but emphasizes that 
the total length of the unoccupied region is L — Nl^ where A'' 
vertical thick bars denote the borders of the areas occupied 
by those excitons. The ends of the SWNT are assumed to be 
a cyclic boundary. 



are eventually reduced to A excitons (as an expectation 
value) through such an EEA process. We make the fol- 
lowing two assumptions: (1) Once two excitons intersect 
in a SWNT, EEA occurs instantaneously with a prob- 
ability of one, and (2) the positions where excitons arc 
created in SWNTs by optical excitations are random. 

Figure Hfb) schematically shows a situation in which 
N excitons exist in a SWNT of length of L. Each exci- 
ton is considered to "occupy" a characteristic length 
in a SWNT. In a static limit, should simply be the 
exciton size. On the other hand, when exciton diffusion 
is present, /x is assumed to be the average distance trav- 
eled by an exciton during the spontaneous decay lifetime 
Ttot and is cx V-Drtot, where D (cm^/s) is the exciton 
diffusion constant. Namely, it is assumed that is de- 
termined by the diffusion length, and any two excitons 
created within ^x undergo EEA. 

First, a 1-D space of length L is considered where N 



segments (or excitons) of length are randomly present 
without overlapping each other, as shown in Fig. [H^b). 
The probability of a new segment (of length l^) to enter 
the system without overlapping any of the N existing 
segments, p{N), is given by the product of the following 
two probabilities 



(1) 



which is the probability for the center of the new seg- 
ment to land on an unoccupied area [hatched regions in 
Fig.Hb)], and 



P2 



L-NL 



N 



(2) 



which is the probability that the occupying length (Zx) 
of the new segment (whose center has landed on an un- 
occupied area) does not interfere with any of the N ex- 
isting segments [vertical thick bars in the lower part of 
Fig. IHb)]. Hence, p{N) is written as 



1 - 



I. 



L-NL 



N 



(3) 



The expectation value of the increment of N due to the 
introduction of a new segment into the system, (AiV)^ 
(0 < (AiV)jY < 1), depends on the type of two-particle 
annihilation considered. For the general "ex+ex X ex" 
{AN)j^ is expressed as 



case. 



{AN)^^p{N)-{l-X)il-piN)). 



(4) 



In the following, the derivation for the case of A = 1 is 
shown as a specific example, since this case is considered 
to represent our experimental situation. The final results 
will be presented for both the A = 1 and A = cases. 



B. Solutions in limiting cases 

1. Steady-state limit 

Here, we consider the steady-state limit, where the 
number of excitons N in the system is steady, correspond- 
ing to CW excitation. In order to derive the relationship 
between /pump and IpL, we consider the relationship be- 
tween Gin and Gout [see Fig. lUJa)], where both Gin and 
Gout are rates, having units of s~^. When N is sufficiently 
small and EEA is neghgible. Gin = Gout = 7totA^- When 
EEA is present, however, this relationship becomes 
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7totiV 



In the case of A = 1, Eq. (HI and Eq. ^ lead to 

7totA^ 



Gi, 



N 



(5) 



(6) 
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We introduce two dimensionless variables C = Nl^/L and 
i^i = GinZx/7totij the former of which is a dimensionless 
exciton population (0 < C < !)• Using these variables, 
Eq. ^ can then be rewritten in dimensionless form: 



c 



(i-C){i-^(i-C)"'} 



(7) 



Expanding the second factor in the denominator of 
Eq. ([7|) and eliminating higher-order terms of l^/ L leads 
to an equation of only ^ and ip, expressed as 



c 



(i-C)E^(t^)'' 



(8) 



Summing up to k = 6 in Eq. ([5]) is sufficient to repro- 
duce Eq. d?]) for l^/L < 0.1. Finally, noting that the 
denominator of Eq. ([8]) is equal to the Taylor expansion 
of an exponential function, the solution for the steady- 
state limit for A = 1 {ex + ex — > ex) is expressed as 



c 



1-C 



exp 



c 



1-C 



(9) 



Here, since /pL oc N and /„ 



oc Gin, /pL and /„ 



proportional to C and tp, respectively, i.e., /pL = ciC and 
-fpump = where ci and C2 are real constants. 

On the other hand, the solution for the case of A = 
(ex + ex ^ 0) is derived through a similar procedure, 
yielding 



c 



2(1-C)exp 



-C 



1 



(10) 



Equations ([9]) and pO|) . valid for CW PL experiments, 
are implicit equations relating the PL intensity (/pl) and 
the pump intensity (/pump) in terms of their respective 
dimensionless variables C and V'- These equations contain 
no fitting parameters except the two linear scaling factors 
ci and C2 and simply become equivalent (C — ip) in the 
low density limit (^ ^ 0). 



2. Instantaneous limit 

Here, we consider the instantaneous limit in which cre- 
ation of all the excitons by an infinitesimally short optical 
pulse and their internal relaxation to En level occur in- 
stantaneously at t = 0, before diffusion-limited EEA and 
spontaneous decay processes begin to occur subsequently. 
This limit represents a situation where the duration of 
the optical pulse and the time required for intraband re- 
laxation are much shorter than the spontaneous decay 
time Ttot, as in the case of the present experimental situ- 
ation. The pump intensity /pump in this limit is directly 
proportional to A'o, the number of excitons or e-h pairs 



created at i = 0, while the PL intensity /pL is propor- 
tional to the number of excitons N that survived EEA. 
The relationship between iVg and N is described by the 
differential equation 



dN 



(AN) 



(11) 



N 



where {AN)j^ is given by Eq. Q. When A = 1 is as- 
sumed, Eq. (fTTj) is rewritten as 



dNo 
dN 



(1 



Nl- 
L 



N 



(12) 



As in the previous case (steady-state limit), we introduce 
dimensionless variables ( = Nl^/L and = NqI^/L. 
Again, since /pL cx; N and /pump oc Nq, /pl = CiC and 
/pump = C2'0; where ci and C2 are constants. With ^ and 
ip, Eq. (|12[) can be written in a dimensionless form: 

_ 1 



(i-C){i-^(i-C)"'} 



(13) 



Furthermore, similar to the steady-state case, an expan- 
sion is performed on the second factor in the denominator 
of Eq. together with elimination of the higher-order 
terms of l^/L, resulting in a differential equation of only 
C and ip- 



dtp 



1 



1-C 



exp 



c 



1-C 



(14) 



Finally, by integrating Eq. from to C, we obtain 
the solution for A = 1 (ex + ex ^ ex) as 



Ei(i; 



where 



Ei(j;) 



y 



-dy 



(15) 



(16) 



is the exponential integral. 

Similarly, the solution for A = {ex -I- ex — > 0) in the 
instantaneous limit is 



2(1-C)exp(-^) 



- 1 



-dC 



(17) 



where the integral has to be solved numerically. Those 
dimensionless equations (|15p and (jl7p become equivalent 
(C = V') in the low density limit (C ^ 0). 



IV. COMPARISONS AND DISCUSSION 

A. Comparison of model and experiment 

To compare with the experimental data shown in 
Fig. [31 we use Eq. which is for the instantaneous 



limit with A = 1. The choice of this equation is because 
of the short duration of the optical pulses used (~250 fs), 
the very fast (^40 fs^^) and efficient^li^S internal decay 
of excitons created at E22 levels to the lowest En level, 
and the much longer spontaneous decay time from the 
Ell level to G. S. (10-100 ps)Md^ The solid curves 
shown in Fig. [3] are fitting curves using Eq. (|15p , indicat- 
ing that the model agrees well with the experimentally 
observed PL saturation curves. 

Table H] is a summary of the fitting analysis performed 
on the data in Fig. [3] using Eq. The first two 

columns show the optimum values of the scaling fac- 
tors Cl and C2~^. These are thought to be propor- 
tional to the oscillator strength for the PL emission 
from Ell levels and the oscillator strength for the op- 
tical absorption around E22 levels, respectively, as ex- 
pected from JpL oc ciTV and A^o ^ C2~^/pump- The right 
three columns (C, N, and R) are values at 1.02 x 10""^* 
photons/cm^ (see Fig. The values of TV were ob- 
tained through the relationship TV = C/^x, where was 
assumed to be 45 nm (one half of the average exciton 
excursion range defined in Ref. [sij). The exciton den- 
sity in the highly saturated regime ('^ 1 x 10^ cm~^) is 
still more than one order magnitude smaller than the ex- 
pected Mott density in SWNTs (~ 7 x 10^ cm~^, assum- 
ing an exciton size of ~ 1.5 nm^), as has been discussed 
in detail in Ref. [l^- R in the right-most column denotes 
the ratio between the number of as-created e-h pairs and 
the number of excitons that survived the EEA process 
until their spontaneous decay to G. S., or Nq/N. The 
values of R show that approximately 90% of the initially 
created e-h pairs/excitons decay non-radiatively through 
the EEA path when the E22 levels are resonantly excited 
at a fluence of 1.02 x lO^"' photons/cm^. 

As expected, the values of ci are similar within the 
same chirality type, regardless of the excitation wave- 
length. However, for the particular case of (6,5) with 
570 nm excitation, where the excitation wavelength ex- 
actly coincided with the E22 resonance peak, the value 
of Cl (~ 3 X 10^) is appreciably smaller than those at the 
other excitation wavelengths (~ 4 x 10^). Such a smaller 
Cl value implies that more excitons decayed to the G. S. 
through non-radiative paths compared to the other cases. 
This may be due to an emergence of stronger nonlinear 
processes that have not been taken into account in the 
assumptions, such as the breakdown of the A = 1 as- 
sumption and/or the appearance of more-than-two-hody 
annihilation processes because of the very high initial e-h 
pair density, achieved with the strong E22 resonance in 
this particular case. 

The values of C2~^ also show the expected tendency 
toward higher (lower) values when the excitation wave- 
length is closer to (further from) the E22 resonance wave- 
length. The slightly higher C2^^ value for the case of (7,5) 
with 570 nm excitation is again attributed to its vicinity 
to the E22 phonon sideband at ^585 nm^ 
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TABLE I: Summary of the optimal parameters to used 
fit the experimental data, obtained from the analysis with 
Eq. (fT^ . (, N, and R are values at the fluence of 1.02 x 10^"' 
photons/cm^, corresponding to the largest fluence in Fig. [3] 
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FIG. 5; (color online). Comparison between the saturation 
behavior predicted by Eq. (fT^ (black sohd curve) and the 
result obtained by a Monte Carlo simulation of the EEA pro- 
cess. The abscissa denotes the density of initially created 
excitons while the ordinate corresponds to the resultant den- 
sity. The experimental range of No in Fig. [3] is No < 2 x 10^ 
cm~^. 



B. Comparison with Monte Carlo calculation 

To further confirm the validity of the model, we per- 
formed a computational simulation based on the Monte 
Carlo method. At the beginning of the simulation, a ran- 
dom distribution of iVo excitons along a line is created at 
time t = 0, corresponding to the instantaneous limit. 
For t > 0, each exciton undergoes a random movement 
in each computational step dt by the distance given by 
the probability distribution of iV(0, Iq), where N{xq, cr^) 
denotes the normal distribution centered at a;o with vari- 
ance (T and Iq = V Ddt. Upon intersection of any two 
excitons, the EEA of {ex -\- ex ex} takes place. In ad- 
dition, each exciton is eliminated from the system with a 
probability of 7tot'^^ (= Ttot~^dt), corresponding to pos- 
sible spontaneous decay during each computational step. 
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The total density of excitons that have decayed sponta- 
neously (N) can thus be calculated for each value of the 
initial density Nq. 

Figure [5] shows a comparison between the model for 
the instantaneous limit [Eq. (jisp ] with = 45 nm and 
the result of the Monte Carlo simulation, plotting the re- 
lationship of A'o and N. The simulation was performed 
with D — 0.42 cm^/s and Ttot = 100 ps,^'^ which resulted 
in an average exciton displacement of 45 nm based on 
the same simulation performed without EEA. The com- 
parison shows satisfactory agreement, indicating that the 
simple analytical form of Eq. well describes the phe- 
nomenon of diffusion-limited EEA, hence validating our 
model. 



Comparison with conventional rate-equation 
solution 



Conventionally, a rate equation of the form 

= Gin(i) - JtotNit) - jEEAN{tf (18) 

at 

has been used to explain EEA processes observed for vari- 
ous materials from 1-D to S-D.^^.J^^iSiSaii^^e j.-^^ ^^^^^ 

on the right hand side, from left to right, represent the 
rate of excitons entering the system (or the En level in 
the present case) , the rate of excitons spontaneously de- 
caying from the system to the G. S. (at 7tot)j and the rate 
of excitons leaving the system because of EEA (at 7eea), 
respectively. Here, the probability of finding two excitons 
at the same position (or the EEA rate) is assumed to be 
proportional to iV^, and hence, exciton diffusion as well 
as the finite size (i.e., length) that each 1-D exciton oc- 
cupies are not taken into account in Eq. (jlSp . 

Equation (fT5|) is solved for the pulse-wise creation of 
iVp excitons at i = [i.e., Gin{t) = S{Na)] in order to be 
compared with the proposed model for the instantaneous 
limit [Eq. (|15|) ] as well as with the experimental results. 
With this initial condition, Eq. ^TE\\ is readily solved to 
be 



N{t) 



+ r exp (7toti)-r 



7EEA 

7tot 



(19) 



The total number of photons emitted from the sample, 
/pL, is obtained by integrating Eq. (fT9| from t = to oo 
as 



'PL = ?7 / "ftotN{t)dt 

'0 



r 



In n - 1 



1 



(20) 

where ry = jr/jtot as before. Equation ([20|) describes 
the relationship between /pL and A'o in terms of two ad- 
justable components rj/T and NqT. 

The dashed lines drawn in Fig. [3] are the optimum fits 
by Eq. to the experimental data. While Eq. (1^ 
reproduces the behavior well for the regime of weaker PL 



saturation with non-resonant excitations, the deviation 
becomes clearer for the regime of stronger PL saturation 
with resonant excitations. Most importantly, however, 
the density predicted by Eq. (|20p has no upper limit as 
recognized by the steadily increasing dashed lines shown 
in Fig. [3] [or by the logarithmic form of Eq. ([20]) ] , while 
Eq. p5p correctly demonstrates the existence of an upper 
limit, as observed experimentally and in the Monte Carlo 
simulation. 

Another important difference between Eq. (|15p and 
Eq. (|20|) is described as follows: The solution of the con- 
ventional rate equation [Eq. (|20p ] contains two indepen- 
dent quantities: rj/T and NqT. The former can be used 
to scale /pL- However, since the parameter A^oF simul- 
taneously scales /pump (fx A'o) and changes the shape 
of the saturation curve, one essentially cannot estimate 
the density of excitons from the fitting analysis based on 
Eq. pp]) . On the other hand, since the shape of Eq. (fT5)) 
has been uniquely determined by its dimensionless rep- 
resentation, the only thing one can do is to linearly and 
independently change the two scaling constants ci and C2 
for /pL and /pump, respectively, yielding sets of C and tp 
values along the scaled (or fitted) /pL vs. /pump curves. 



V. SUMMARY 

We have investigated photoemission properties of high- 
density 1-D excitons in single-walled carbon nanotubes. 
As the excitation intensity increases, all photolumines- 
cence emission peaks arising from different chiralitics 
showed clear saturation in intensity. Each peak exhibited 
a saturation value that was independent of the excitation 
wavelength, indicating that there is an upper limit on the 
exciton density for each nanotube species. We interpret 
these results in terms of diffusion-limited exciton-exciton 
annihilation processes through which high-density exci- 
tons decay non-radiatively. 

To quantitatively understand the saturation behavior 
observed in the experiment, we have developed a theoret- 
ical model, taking into account the generation, diffusion- 
limited exciton-exciton annihilation, and spontaneous de- 
cays of 1-D excitons. The saturation curve predicted by 
the model under appropriate approximations reproduced 
the experimental saturation curves well, and the fitting 
analysis allowed us to estimate the density of excitons for 
a given diffusion constant. We also compared our results 
with Monte Carlo calculations, confirming the validity 
of our model. Additionally, we examined the saturation 
behavior predicted by the solution of the conventional 
exciton-exciton annihilation rate equation that does not 
take into account diffusion and showed that the solu- 
tion qualitatively failed to fit the experimental data at 
high exciton densities, showing the inappropriateness of 
its use for excitons in carbon nanotubes. The approach 
presented in this paper should have wide applicability 
for predicting the intensities of photoluminescence from 
single-walled carbon nanotubes under various excitation 
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conditions. 

An important conclusion drawn from the current study 
is that a large density of electron-hole pairs is difficult to 
achieve in single-walled carbon nanotubes. This imposes 
a serious challenge for making active optoelectronic de- 
vices based on semiconducting single- walled carbon nan- 
otubes — 1-D semiconductor lasers, in particular, which 
would require densities comparable to, or higher than, 
the Mott density. The existence of an upper limit on 
the density of excitons would also prevent fundamental 
studies of bosonic characters of 1-D excitons expected 
at quantum-degenerate densities. Hence, novel methods 
are needed for minimizing non-radiative decay processes 
such as exciton-exciton annihilation in single-walled car- 
bon nanotubes, which would require ways to control the 
dynamic parameters of excitons such as the spontaneous 
decay rate and diffusion constant through, e.g., varying 



the temperature as well as applying an external magnetic 
field41 
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